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Abstract 



With the global color symmetry model being extended to finite chemical po- 
^ ■ tential, we study the density dependence of the local and nonlocal scalar quark 

condensates in nuclear matter. The calculated results indicate that the quark con- 
densates increase smoothly with the increasing of nuclear matter density before the 
critical value (about 12po) is reached. It also manifests that the chiral symmetry is 
restored suddenly as the density of nuclear matter reaches its critical value. Mean- 
while, the nonlocal quark condensate in nuclear matter changes nonmonotonously 
against the space-time distance among the quarks. 
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1 Introduction 



It is well known that quark condensates are the essential characteristics in the phase 
transition from the quark gluon plasma to hadrons of the early universe evolution. Mean- 
while the condensates feature the nonperturbative structure of the QCD vacuum, and 
hence play essential roles in describing hadron structure, and further the properties of 
nuclear matter and finite nuclei [TJ El HD lU El E]- In the aspect of describing hadron 
structure in QCD, not only the local scalar quark condensate <: qq :>, but also the 
nonlocal condensates <: q(x)q(0) :> are required. Especially, the later is more im- 
portant to figure the nonpoint particle property of hadrons jlj [7J El El QUI- Then a 
large body of works on the quark condensates in both free space and hot and dense 
nuclear matter have existed (for a recent review, see for example Ref.|llj). Many of 
the approaches give a descent feature of the local scalar condensate in nuclear matter 
at low energy, which is believed to be the manner for the chiral symmetry to be re- 
stored. However an "upturn" emerges at higher density in the linear Walecka model[12j, 
Dirac-Brueckner method JH] and Schwinger-Dyson formalism |14j. Recently, the Dyson- 
Schwinger equation (DSE) method of QCD ^3^21 and the global color symmetry model 
(GCM) of QCDjT71 EB1 HH1 ED] have also been employed to evaluate the condensates in 
vacuum [2TJ E21 EH| and those at finite temperature [I] and finite density [2H ES]- The DSE 
method calculations |24[ 123] show that the quark condensate increases with respect to the 
increase of chemical potential (or the nuclear matter density) before a critical value is 
reached. It means that one has not yet reached a common realization on how the quark 
condensate changes against nuclear matter density qualitatively. Since the DSE approach 
and the GCM have been shown to be quite successful in describing hadron properties 
both in free space and in nuclear matter [2H E3 EH E3 EH EH EH], with the global color 
symmetry model being extended to that at finite nucleon density p (as usual, it is achieved 
by involving a finite chemical potential /i), we will study the dependence of the local and 
nonlocal scalar quark condensates on nuclear matter density again in this paper. 

The paper is organized as follows. In Sec. II we describe the formalism of the GCM at 
finite chemical potential /i and the relation between the chemical potential and the baryon 
density in nuclear matter. In Sec. Ill we represent the calculation and the obtained 
results of the local and non-local scalar quark condensates as functions of the nuclear 
matter density. In Section IV, a brief summary and some remarks are given. 
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2 Formalism 



The starting point of the global color symmetry model (GCM) is the action in Euclidean 
matric[T7j 

,2 



S = Id xd y 



q(x) (7 • d + m) 5{x - y)q(y) - %-jl(x)D^ u (x - y)j"(y) 



where j° = q{x)^-^f ll q{x) is the quark current, D fMU (x — y) is an effective two-point gluon 
propagator, m is the current quark mass, and g is the quark-gluon coupling constant. 



Taking the effective gluon propagator to be diagonal, i.e., D 



X 



V) 



8 llv D(x 



y) and 



applying the Fierz reordering to the quark fields, one can rewrite the current-current term 



as 
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d*xd*yj a Jx)D(x - y)f v {y) = / d*xd?yj°(x, y)D(x - y)j°(y, x) 



2 J v- 2 

where j (x, y) = q(x)A e q(y) with A e being the direct products of Lorentz, flavor and color 
matrices which produce the scalar, vector, pseudoscalar and axial vector terms. More 
concretely, A 8 = \K a ® C b ® F c with {K a } = {I, z 7s , J= 7m , ^j 5 }, {C b } = {f J, ^=A}, 



{F c } = {^=1, ^} where {X a }(a = 1, 2, • • • , 8), {f }(i = 1, 2, 3) are the generators of the 
groups SUc(3) and SUp(2), respectively. It is obvious that such a j (x,y) is a bilocal 
current. With two flavors u and d of quarks being taken into account, each A e is either 
isoscalar or isovector. The color matrices involved in the Fierz transformation contain 
color-singlet and color-octet terms. Taking the bosonization procedure one can transfer 
the bilocal quark current structure into auxiliary Bose fields carrying the quantum number 
6. The action of the GCM in free space (i.e., at the chemical potential jj, — 0) for the 
zero-mass quark can then be rewritten ^7] in the Euclidean space as 



S(B) = / d 4 xd 4 yq(x) 7 • d5{x - y) + A e B e (x, y) q{y) + 



4 4 B e (x,y)B e (y,x) 
V 2gW(x-y) ' 

where B 9 (x,y) is the bilocal Bose field. 

To extend the GCM to finite nuclear matter density (with finite chemical potential 
iu), one should, in view of the statistical mechanics, take the partition function of the 
canonical (quark) ensemble into that of the grand ensemble with quarks and hadrons that 
are the solitons collecting quarks [THj. The quark field should be transformed under a 
constraint on the baryon number through the chemical potential \x 



q(x) 



q (x) 



q(x) 



(1) 



After some derivation, we have the action of the GCM in nuclear matter 

S(B 8 ,^) = fd 4 xd 4 yq'(x)\(-fd-fx l4 )5(x-y) + e> 1Xi A e B e (x,y)e-^ 4 } q'(y) 



+ J d A xd A y 



B%x,y)B e (y,x) 
2gW(x - y) 
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and the generating functional is given as 

Z[/i, rf,rj\=J DqDqDB e ^- s(B9 > ti+J^im+lv)] (2) 
where rj and r\ are the quark sources. After integrating the quark fields, we obtain 

S(B e , fj) = -Trln \( 1 -d-fi^)5(x-y)+e^A B e e-^] + [ d A xd A y B ^ V ] B ^ . (3) 

J J 2q 2 D(x — v) 

writte 

'X + y 



Generally, the bilocal field B (x,y) can be written asj 

B\x,y) = B„ 9 (^) + Er!(x,!/)4^) , (4) 



where B^x^y) =Bq(x — y) — B e (x — y) is the vacuum configuration of the bilocal field. 
S^rf (x, y)4>i {^2 L ) correspond to the fields which can be interpreted as effective meson 
fields. In the lowest order approximation with only the Goldstone bosons being taken into 
account, the 0f includes n and a mesons. The corresponding width of the fluctuations 
are approximately the same as the vacuum configuration 20], i.e., Tf = Bq. Since the 
bilocal field arises from the bilocal current of quarks, the internal q-q structure of the 
mesons can be well described in the GCM. The vacuum configuration can be determined 
by the saddle-point condition = 0. Then an equation of the translation invariant 
quark self-energy E(g, /x) is obtained as a truncated Dyson-Schwinger equation 

f d A q t a 1 t a 

E (P^) = / 77T^l9 2 D(p- q)—lu- -=7 sl^i (5) 

J (27r) 4 2 ry ■ q - /ry 4 + T,(q, fi) 2 

Considering the fact that the inclusion of the chemical potential breaks the 0(4) symmetry 
in the four-dimensioanl space, one should rewrite the decomposition of the self-energy E 

as 

V) = i[ A (P, AO - 1]7 * P + *[C0> ^) - 1]74(P4 + ¥) + B (p, /1) , (6) 
where B(p,/j,) is the counterpart of the vacuum configuration in the momentum space, 



i.e., 



B{P^) = J -±- 4 B d (z)e^dz. 



(2nY 

It requires then that the vacuum configuration of a and 7? should satisfy the restriction 
a 2 + tt 2 = 1. With q u = (q,q4 + ifi) being introduced, and combining Eqs. (5) and (6), 
one can obtain the equations to determine the A(p), B{p) and C(p) as follows 

\A(~\ 11 -B 8 f d ^ 2n/ \ A(q)q-p 

=3]^ Dip - q) A 2 (qW + C 2 (q)q 2 + B 2 (q) ' ^ 

r ^/~x ,1-2 8 f A 2 C{q)q 4 p 4 

[C{P) - 1]P *=3Jj2^ 9 D{P ~ q) A>(q)? + C*(q)ft + B*(q) ' (76) 
16 f d*q 2 B(q) 

B{P) = YJ(2^ 9D{P - q) A 2 (q)q- + C 2 m + B 2 (q) ' ^ 



4 



Basing on the solution of the Dyson-Schwinger equation or the gap equations, one can 
determine the bilocal field, and fix further the GCM action S and the quark propagator 
G with 

G- 1 = iA{p)^ ■ p + iC{p) lA p A + A e B\p) . (8) 

Thereafter we can evaluate the quark condensates from the above saddle-point expansion. 
It is straightforward that, under the mean-field approximation, for any quark operator 

one can get the condensate 

<: O n :> = (-l) n E Tr [(-) Pr S 1 r 2L • • ■r{l(G)i 1J -P(i)(G) iai p (2 ) • • • (G)i n jP(n)] , 
P 

where G is the quark propagator and P denotes the permutation of the n indices. We 
obtain then the local scalar quark condensate as 

<: qq :>=<: q(0)q(0) :>= -Tr^ c G{x,x) 

_ _o,/2 f ^f) (q) 

J {2TiYfA\p 2 )+piC 2 {p 2 )+B\p 2 ) ' y > 

Meanwhile, the nonlocal quark condensate is given as 
<: q{x)q(0) :>= -Tr jC [G(x, 0)] 

= _o./o f Jpx (in) 

J (27v) 4 fA 2 (p 2 )+plC 2 (p 2 ) + B 2 (p 2 ) y ' 

It is apparent that, with the solutions of the gap equations (7a)-(7c) being taken as the 
input for Eqs. (9) and (10), the variation of the quark condensates against the chemical 
potential can be evaluated. In practical calculation, since the knowledge about the exact 
behavior of g 2 and D(p — q) in low energy region is still lacking, one has to take some 
approximations or phenomenological form to solve the gap equations. For simplicity, we 
adopt the infrared dominant form j3H IT7] 

g 2 D(p-q) = ^5(p-q) 1 (11) 

where rj is an energy-scale parameter and can be fixed by experimental data of hadrons. 
Although this form does not include the contribution from the ultraviolet energy region, 
it maintains the main property of QCD in low energy region. With Eqs. (7a)- (7c) and 
(11), one has the Nanbu-Goldstone solution 

2 

A{p) = C{p) = 2, B{p) = {rf - 4P 2 ) 1 / 2 , for Re(p 2 ) < \, (12a) 



A(p) = C(p) = ~ 



V v 



2 V 2 \ 1/2 ' 



4 

2 

B{p) = 0, for Re(p 2 ) > -j, (12b) 
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which describes the phase where chiral symmetry is spontaneously broken and the dressed- 
quarks are confined. Meanwhile one has also the Wigner solution 



MP) = c(p) = - 



\ + 2 4 

p 



1/2 



B(p) = 0, 



(13) 



which characterizes a phase in which chiral symmetry is not broken and the dressed-quarks 
are not confined. 

To fix the energy-scale parameter rj so that the above mentioned calculations can be 
accomplished, we take the way described in the original work of the GCM[T7|. Extending 
the relation between the GCM and the bag models given in Ref . [T7j to that in nuclear 
matter, we have the relation among the nucleon radius R, nucleon mass M and the bag 
constant B of a nucleon as 



R 



(14) 



An / AT]-/? \ X /4 

where a = 3(u — ji) — Z = 6.12 — 3/x — Z , with Z being a parameter to include the 
corrections of the motion of center-of-mass, zero-point energy, and other effects of the 
three quarks in a nucleon. Meanwhile, in advantage of the Friedberg-Lee nontopological 
soliton ansatz|32], the bag constant comes from the difference of the action for hadrons 
and that for the vacuum. We have then 

B = S'(hadron) — S'(vacuum). 

In the GCM, the interaction is transferred by the Goldstone bosons n and the scalar 
meson a. From the restriction n 2 + a 2 = 1, one can take a simple approximation n = 0, 
a = 1 for the vacuum. Referring the configuration of the mesons in a hadron as (a, 7?), 
we have 

B = S(a, If)- 5(1,0). 

In principle, the configuration of a and n should be determined by solving the coupled 
equations of motion of the mesons and the quarks. As an approximation jT7j, they can be 
taken as a = 0, n = 0. After some derivation we obtain 



B = 5(0,0,) -5(1,0) =12 J-0j-A^ 



A 2 (p 2 )f+C 2 (p 2 )Pl+B 2 (p 2 ) 
A 2 (p 2 )p 2 + C 2 (p 2 )Pl 
B\p 2 ) \ 



A 2 (p 2 )p 2 + C 2 (p 2 )Pt+B 2 (p 2 ) 



(16) 
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In order to investigate the dependence of the quark condensates on the nuclear matter 
density p explicitly, we must transfer the above /^-dependence to the p-dependence. We 
should then determine the relation between the nuclear matter density p and the chemical 
potential p in the GCM. It has been known that the baryon number in nuclear matter 
can be related with the generating functional of the system [T£j ESI IS3J as 

n = ^-hx Z(ji, B e {x-y))^ ^Tr In CT V, B\x - y)) , (17) 

where G~ 1 (p, B e (x — y)) is the inverse of the quark propagator in the medium. Taking 
some algorithm, we obtain finally 



dp J (27r) 4 7 * p 2 
The baryon number density can thus be given as 

Pn fd*x (2^8(1 J P p 2 ' 1 j 

3 Numerical Result and Discussion 

Taking the energy scale by fitting the properties of ir and p mesons rj = 1.37 GeVj21] 
and calibrating the nucleon mass M = 939 MeV in free space in the GCM, we obtain 
the center-of-mass motion correction parameter Z Q = 3.707 and the nucleon radius Rq = 
0.69 fm. With the fitted parameters and the above formulae, we get at first, in the 
Nanbu-Goldstone phase, the quark condensate in free space <: qq :>o= —(148 MeV) 3 . 
Meanwhile, in the Wigner phase, the quark condensate vanishes, i.e., <: qq :>o= 0. It 
is evident that the presently obtained <: qq :> in the Nanbu-Goldstone phase is very 
close to the results given with a fully dressed gluon propagator (— <: qq :> 1 / 3 can be 
150-180 MeV)j2U I2U I2Z| even though it is smaller than those given in QCD sum rules 
and some other approaches. From Eqs.(7) and (10) one can know that the value of the 
scalar quark condensate in vacuum depends on the integration interval determined by the 
parameter rj. With the increase of rj, the absolute value of the condensate will increase. 
In the present calculation, the r\ is fixed by fitting the nucleon properties consistently but 
not freely. What we are now interested in is the changing characteristic of the condensate 
in nuclear matter against the nuclear matter density, which can be identified as the ratio 
of the condensates in nuclear matter to that in vacuum. Since the condensate in nuclear 
matter depends also on the parameter rj with the same relation, the smaller absolute value 
in vacuum will not affect the changing feature. 

Recalling the property of the quark confinement phase and deconfinement phase, one 
can know easily that, in the confinement phase, the motion of a quark is restricted in a 
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hadron. In a phenomenological word, the motion of the quark is limited in a bag. The 
bag constant provides the constraint to limit the motion of the quark. It is apparent that, 
if the bag constant B becomes zero, the quark can move freely. In view of the variation 
of the bag constant, B = can then be taken as the critical point of the confinement 
phase and the deconfinement phase. Consequently, the gap equations (7a)-(7c) take the 
Nanbu-Goldstone solution or Wigner solution. Accomplishing the integral in Eq. (16) 
with various values of the chemical potential p, we get the chemical potential dependence 
of the bag constant. With the relation given in Eq. (18), we obtain the nuclear matter 
density dependence of the bag constant. The result is illustrated in Fig. 1. The figure 
shows evidently that the bag constant decreases with respect to the increasing of the 
nuclear matter density, and a critical nuclear matter density p ~ 11.8po (corresponding to 
a critical chemical potential \i c ~ 0.316 GeV), where po is the normal density of nuclear 
matter, exists for the quark deconfinement to happen. 

By varying the chemical potential p and accomplishing the calculation in Eqs (9), 
(10) and (18) at several space-time distances x, we obtain the variation behavior of the 
ratios of local and non-local scalar quark condensates in nuclear matter to the local 
condensate in vacuum. The results with several fixed space-time distances are illustrated 
in Fig. 2 [simply denoted as R 2 (p)]. The figure shows clearly that, in the Nanbu-Goldstone 
phase, the local scalar quark condensate and the nonlocal condensates with each fixed 
space-time distance between the quarks increase smoothly with respect to the increase of 
nuclear matter density. As the density of nuclear matter reaches its critical value p c , the 
condensates vanish suddenly. In the Wigner phase, the condensates maintain zero since 
B(p) = 0. This result provides a clue that the chiral symmetry is broken more strongly as 
the nuclear matter density increases gradually. When the nuclear matter density reaches 
the critical value, the chiral symmetry is restored suddenly. Such a varying behavior 
is consistent with the result of the DSE calculations j2H 122] and the sudden changing 
characteristic is similar to that given in the composite-operator formalism of QCD 35J. 
However, it is evident that such a changing behavior of the scalar quark condensates in 
nuclear matter is different from most of the previous results. It has been conventionally 
known that the nucleon mass is a monotonous function of the scalar quark condensate, and 
the mass of a nucleon in nuclear matter decreases with respect to the increase of nuclear 
matter density. Then the absolute value of the scalar quark condensate should decrease as 
the nuclear matter density increases. In fact, the condensate depends also on the number 
of nucleons in nuclear matter which increases rapidly against the increase of nuclear matter 
density. Combining these two aspects together, one can know that it is reasonable for 
the condensates to increase with respect to the increase of nuclear matter density. On 
the other hand, according to the well known Gell-Mann— Oakes— Renner relation[36j, the 
mass of a pion in nuclear matter will increase with respect to the increase of nuclear matter 
density. Consequently, the threshold of the pion production in ultrarelativistic heavy ion 
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collisions will increase. One of the reason for the most recently observed relatively small 
amount of pion production to proton production in the RHIC experiment {3 7] may be such 
an increase of the pion production threshold. 

Looking over Fig. 2, one may also easily realize that, at a small space-time distance, 
the condensate ratios at a definite density of nuclear matter are larger than that of the 
local condensate. However, at large space-time distance, they are smaller. It indicates 
that the quark condensate in nuclear matter does not change monotonously with respect 
to the space-time distance. To show this point more explicitly, we display the changing 
feature of the ratio of scalar quark condensate in nuclear matter to the corresponding 
values in free space against the space-time distance in Fig. 3 [denoted as Gr 2 (£ 2 )]- The 
figure shows evidently that the condensate in free space (p = 0) decreases gradually 
with the increase of the space-time distance among the quarks definitely. However, the 
condensate in nuclear matter(p ^ 0) increases if the space-time distance is very small. As 
the distance is about 0.2-0.24 fm (about one fourth of the nucleon radius), the condensate 
ratio reaches its maximum if the nuclear matter density changes from 3po to llpo> an d such 
a distance increases with the increase of nuclear matter density. If the space-timer distance 
gets larger, the condensate ratio decreases with the increasing of space-time distance and 
the decreasing rate increases with respect to the increase of nuclear matter density. It 
manifests that the chiral symmetry can be restored more rapidly with the increase of 
nuclear matter density if the space-time distance among quarks is approximately equal to 
or larger than the one forth of the nucleon radius. Furthermore, such a nonmonotonous 
changing characteristic also indicates that there exists a repulsive interaction at very small 
space-time distance between the quarks. Since increasing the density of nuclear matter 
results generally in a decrease of the distance among quarks, the decrease of condensate 
ratios as the quarks are located farther away beyond the repulsive core is just the other 
manifestation of the fact that the absolute values of the scalar quark condensates increase 
against increasing nuclear matter density. 

4 Summary and Remarks 

In summary we have investigated the nuclear matter density dependence of the local and 
nonlocal scalar quark condensates in nuclear matter in an effective field theory model of 
QCD, namely the global color symmetry model (GCM). The calculated results indicate 
that the quark condensates increase smoothly with the increase of the density of nuclear 
matter and vanish suddenly as the nuclear matter density reaches its critical value (about 
11.8po)- Meanwhile the nonlocal condensates change nonmonotonously as the quarks are 
separated far away from each other. This result manifests that the chiral symmetry is 
broken more seriously as the density of nuclear matter increases and can be restored 
suddenly at the critical point. However, if the quarks are separated farther away from 
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each other, the chiral symmetry restoration process can be enhanced with the increase of 
nuclear matter density. 

Since only a very simple form of the g 2 and D(p — q) is taken into account in the 
present calculation, detailed effects of the running couple constant g 2 (or a), the gluon 
propagator D(p—q) and the other degrees of freedom have not yet been included. Besides, 
the dependence of the quark, meson states on the nuclear matter density (or chemical 
potential) has not yet been taken into account self-consistently. Moreover, even though 
the presently obtained result is coincident with that obtained in the DSE formalism |24[ 12*5] 
and the sudden change feature is consistent with that given in the composite-operator 
formalism of QCD[[J5], the changing behavior in the continuous region is not consistent 
with most of the previous results (e.g., those given in Ref.fll]). Then more sophisticated 
investigation is necessary and under progress. 
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Figure 1: The variation of the bag constant with respect to the nuclear matter density 
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Figure 2: Calculated ratio of the two-quark condensate in nuclear matter to that in free 
space at several fixed space-time distances 
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Figure 3: Calculated ratio of the nonlocal two-quark condensate in nuclear matter to that 
of the local one 
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